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APPROXIMATE ANALYSIS OF THE PERFORMANCE 


OF CHAR-FORMING ABLATORS 
By Robert T. Swann 


SUMMARY* 


Equations governing the performance of charring ablators, subject to a 
diffusion- controlled oxidation mechanism of char removal, are derived. A solu- 
tion is presented for the case of a constant-enthalpy level. The effects of 
various material properties and environmental parameters are discussed. 


INTRODUCTION 


The charring plastics provide a practical method for protecting the inte- 
rior of a reentry vehicle from the aerodynamic heating encountered during 
reentry. Much effort has been directed toward understanding the performance of 
these materials. Analytical models for charring ablators are presented in ref- 
erences 1 and 2 . These references also present estimates of the weight of the 
charring ablator which is required to provide adequate thermal protection as a 
function of heating conditions and material properties. A more rigorous ana- 
lytical model of the charring ablator is presented in reference 3. Results of 
some typical calculations are given in references 3 to 5. 

A primary problem encountered in analyses of the performance of charring 
ablators is the removal of char from the surface. Aerodynamic shear, thermal 
stress, differential pressure, and oxidation have been suggested as mechanisms 
contributing to char removal. The removal of char as a result of pressure 
within the char, caused by confinement of the gaseous products of pyrolysis, is 
analyzed in reference 6. The effects of aerodynamic shear and thermal stress 
are not known in quantitative terms. Results of an extensive experimental 
investigation of the effects of oxidation on char removal and surface heating 
are presented in reference 7- In the test environment employed in the inves- 
tigation (ref. 7), oxidation appeared to be the only mechanism involved in char 
removal. In the light of the results presented in reference 7 * it appears that 
oxidation effects should be incorporated in the analysis of charring ablators. 

The differential equation governing the performance of a charring ablator 
subject to a diffusion controlled mechanism of oxidation has been derived 
herein and is solved for a particular case. The effects of material properties 
and environmental factors on performance are discussed. 
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concentration of oxygen 

concentration of oxygen external to boundary layer 
volatile fraction in uncharred material 
heat of combustion per unit weight of char consumed 
effective heat of pyrolysis 
heat of pyrolysis 

total enthalpy 
thermal conductivity 
rate of mass transfer 
rate of char removal 

mass of gaseous products of pyrolysis 
rate of formation of gaseous products of pyrolysis 
Lewis number, %ij/^Sc 
P randtl number 
Schmidt number 

heating of surface resulting from combustion 
hot-wall convective heating rate 

hot^wall-. convective heating rate with no mass transfer 

radiant heating rate 
heat input 
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Q c total convective heat input 

Q r total radiant heat input 

t time 

T ab s olut e t emp eratur e 

T s average surface temperature during reentry 

u boundary- layer velocity parallel to surface 

v boundary- layer velocity normal to surface 

W thermal-protection weight per unit area (weight of material degraded) 

W c char weight 

x boundary- layer coordinate parallel to surface 

y boundary- layer coordinate normal to surface 
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char thickness 
emissivity 

transpiration factor 

weight of char removed per unit weight of oxygen diffusing to surface 

viscosity 

char density 

density of boundary-layer fluid 
Stefan-Boltzmann constant 


a 

Subscripts: 

e external to boundary 

i initial value 

m extremum value 

max maximum value 



0 value in absence of mass transfer 

p pyrolysis 

r radiant 

ref reference value 

s surface 

w at the wall 


ANALYSIS 


An analytical model of the charring ablator is presented in reference 1. 
This model is shown schematically in figure 1. The outer surface is subjected 

to aerodynamic heating. The char pro- 
vides an insulating layer, and the heat 
passing through this layer is absorbed 
by pyrolysis at the interface between the 
char layer and the uncharred material. 

The surface of the char re-radiates a 
significant fraction of the heat input. 
The gases generated by pyrolysis tran- 
spire through the char layer and are 
injected into the boundary layer. The 
gases are heated as they pass through 
the char layer, and this removal of heat 
from the char layer reduces the quantity 
of heat conducted to the pyrolysis inter- 
face. When these gases are injected into 
the boundary layer, the convective heat transfer to the surface is reduced. 

This reduction in convective heat as a result of mass injection is the same 
effect as that obtained with simple subliming ablators. In addition to the 
gases produced by pyrolysis, a residue remains at the interface and adds to the 
thickness of the char layer. This process tends to increase the thickness of 
the char layer. At the same time, char removal may occur as a result of thermal, 
chemical, or mechanical processes, and tends to decrease the thickness of the 
char layer. 

The thermal behavior of charring ablators is very complex and, in general, 
numerical methods are required to solve the governing equations. However, if 
some simplifying assumptions are made, approximate expressions can be obtained 
for the char accumulation and the pyrolysis rate. Such approximate results are 
compared with the numerical results of reference 1, and the agreement is satis- 
factory. The effects of oxidation on char removal and heat input can be 
included in analyses similar to that of reference 1 without any particular dif- 
ficulty. However, equations must be obtained for the rate at which oxidation 
removes char and the rate at which oxidation adds heat to the surface. 


(Reradiation) m p q c + q p m c q 



i 'O' 1 

l 

? 

Char layer 

"p 

? ? ? 

1 


Uncharred 

material 



Front 

surface 


Pyrolysis 

interface 


Figure 1.- Schematic diagram of charring 
ablator. 
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Oxidation 


It is shown in reference 8 that the rate of oxidation of porous carbon in 
air is a diffusion controlled process, even at relatively low temperatures. 
Therefore, the rate of oxidation is controlled by the rate at which oxygen dif- 
fuses from the external stream to the surface. A detailed study of the removal 
of char by oxidation is presented in reference 7* I n that study it is concluded 
that char removal occurs primarily as a result of oxidation. With sufficiently 
high surface temperatures (and hence high reaction rates), the oxygen diffusing 
to the surface reacts with the char rather than with the products of pyrolysis. 
The maximum rate at which oxygen diffuses to the surface is 


m(02) = 


C 

u e iN Le 


He “ H w 

The derivation of equation (l) is outlined in the appendix. 


( 1 ) 


If A is the weight of char removed per unit weight of available oxygen, 
the rate of char removal is 


m n — 


c n ■, _ f ic 
c e N Le ^ - a 


He ~ H w 


H e - H- 


( 2 ) 


■w 


If Ah is the quantity of heat added to the surface per unit weight of char 
consumed, the rate at which the surface is heated by combustion is 


qi H e - H w 


a Ah 


( 3 ) 


The evaluation of A and Ah is a complicated aero- thermal- chemical 
problem, the solution of which is beyond the scope of the present paper. How- 
ever, limiting values can be indicated. The maximum value of A is 0.75> 
which corresponds to the ratio of the molecular weights of carbon and oxygen 
when the product of oxidation is carbon monoxide. For this reaction, the value 
of Ah is approximately 000 Btu/lb. Carbon monoxide is the most probable 
product of a reaction occurring in an oxygen-deficient environment, and the 
formation of it is consistent with the results of reference 7- 

The model of a charring ablator described previously is adequate for numer- 
ical analysis. With certain simplifications, the model can be reduced to a form 
for which analytical solutions can be obtained. First, it is assumed that the 
specific heat of the char layer is negligible. Solutions based on this assump- 
tion are compared with numerical solutions in reference 1 and the approximate 
solution provides an adequate representation of phenomena within the char, 
including char formation. A second assumption is that the effective conductive 
ity of the char layer is a function of temperature as follows: 


5 


ii 




Because the specific heat of the char layer is negligible, the heat transferred 
across the char layer (from the surface to the reaction zone) is equal to 
k(T s - T p ) 

— * It is expected that the thermal conductivity of a porous material 

^c 

such as a char layer will increase with increasing temperature because of the 
greater importance of radiation across the void spaces. The surface 
re-radiates energy at a rate proportional to its fourth power; therefore, rela- 
tively large changes in the heat input are required to obtain significant 
changes in the surface temperature. Therefore, by properly selecting the refer- 
ence temperature T re ^, the error introduced by using equation (k) can be mini- 
mized, The actual error cannot be estimated because no data are presently 
available. 

The model of a charring ablator presented in reference 1 can be readily 
extended to the case in which char removal and surface heating occur as a result 
of oxidation. The unknown quantities are the surface temperature, the rates of 
pyrolysis and char removal, the convective and combustive heating rates, and 
the weight of char accumulated. From these quantities, the weight of material 
degraded can be derived. 

The approach which will be used is as follows. Energy balances at the 
heated surfaee and at the pyrolysis interface are used to eliminate the surface 
temperature. An approximate boundary- layer solution is used to express the con- 
vective heating rate in terms of the convective heating rate with no mass trans- 
fer and the mass-transfer rates. The char- removal rate and the combustive- 
heating rate are related to the convective-heating rate by equations ( 2 ) and ( 3 ). 
An equation for the weight of char accumulated is obtained by noting the com- 
bined effects of char formation and char removal. A differential equation 
relating the weight of char accumulated and the convective heat input with no 
mass transfer is then derived. This equation, together with the relation 
between char weight and amount of material pyrolyzed, is used to relate the 
weight of material degraded and the convective-heat input with no mass transfer. 

A heat balance at the outer surface yields the following equation: 

% + «1 + 1 r = ^ + |£(t s - Tp) 

C (5) 

Convective Combustion Radiant Re- radiation Heat conducted 

heating heating heating from surface 

The heat conducted from the surface is absorbed by pyrolysis and by heating 
the gaseous products of pyrolysis to the surface temperature. The rate of 
pyrolysis is 
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(6) 


M T s - T p) 
w c ^eff 


where 


Ah eff = Ah p + ^(T s - T p ) (7) 

is an effective heat of pyrolysis which includes the heat absorbed by the gases 
as they transpire through the char. Equation (7) is derived in reference 2. 
Equations (4), ( 5 ), and (6) can be combined to eliminate T 


q r + 


qi + q c = ip Aheff 




(8) 


where 


p = ^ 

ae(T re f 2 + T p 2)( T ref + T p) 

and it is assumed that 

T v k « T s 4 


The convective-heating rate can be obtained from the relation 


Q-c 3-c, o 


1 - Ti(m c + ip) 


H, 




l c.o J 


(9) 


where q c Q is the hot-wall convective-heating rate with no mass transfer. 

Equation (9) is obtained by fitting a straight line to the boundary- layer solu- 
tions given in reference 9* It is valid only when the term 



+ m 


*p) 


H e - ^ 


< 


2 

3 


This condition is almost always satisfied when all the heat input is convective. 
A higher order approximation may be required if high radiant-heat inputs are 
also incident on the surface. It is assumed in equation (9) that the char which 
is oxidized and the pyrolysis products which are injected into the boundary 
layer are equally effective in blocking the convective heat input; this assump- 
tion explains the multiplication of m c and mp by the same value of tj. 

Actu ally , the blocking effectiveness depends on the molecular weight of the 
injected gases and probably differs for the products of combustion and pyrol- 
ysis. However, to the present order of approximation this difference may be 
neglected. 
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From equations (2) and ( 9 ), the rate of oxidation of the char is 


Ac 


a / ^c, o 
1 ■+■ aT| \H e — H w 



( 10 ) 


and similarly the heat transfer to the surface (in terms of the heat transfer 
with no mass transfer and the rate of pyrolysis) is 


„ _ (H e - H w )£ c _ 1c, o “ ^(He “ H w ) 

qc = a = 1 + at) 


( 11 ) 


From equations ( 3 ) and (ll), the rate at which the surface is heated as a 
result of combustion is 


!l 




a Ah 
1 + at) 


(12) 


To this point, a sufficient number of relations has been developed to 
express all the unknown variables in terms of the rate of pyrolysis, the weight 
of char accumulated, and the heating rate with no mass transfer. From consider- 
ation of the geometric effects of pyrolysis and char removal, the weight of char 
accumulated is related to the rates of pyrolysis and char removal as follows: 


aw c 

dt 



- m r 


( 13 ) 


which states that the rate of change of char weight is equal to the difference 
between the rate at which char is formed and the rate at which char is removed. 
Equation (l0) can be used to eliminate m c from equation ( 13 ). Thus, 


dW c 

dt 


= m T 


1 - f + aq 
_f(l + aq)_ 


aq. 


c , o 


(1 + aq ) (H e - H w ) 


(14) 


Equations ( 8 ), (ll), (12), and (l4) can be combined to yield the fol- 
lowing equation: 


*c,o a Ah + H e - Hy + q r 
H e - H w Ah eff (l + aq) Ah eff 


1 - f + aq 


dW 


(1 + a q)— - + 


C a lc,Q 


dt 


He “ H w 


( a Ah + H e - H w ) 


W c 

+ 1 + T] , 

J3 Ah eff (l + aT l) - 1 


( 15 ) 


The quantities q c Q , q r , and H e - H w are, in general, known functions 
of time. In an important special case to be discussed later, that is. 
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laboratory evaluation of charring ablators, W c is a function of time whereas 
q. c , o? lr> and - H e “ H w are constants. Equation (15) is a nonlinear first- 

order ordinary differential equation which can be solved by numerical methods. 
When the weight of char has been determined, equation (lk) can be integrated to 
determine the mass of the gaseous products formed by pyrolysis. The weight of 
material degraded to yield gaseous products and char is 

w = ^ + w c ,i (16) 


where W c ^i is the char initial weight. 

The equation governing the performance of a charring ablator (eq. (15)) 
cannot be integrated in general terms; however, certain significant information 
can be deduced from the equa- 
tion without obtaining a solu- 
tion. In addition, a solution 
can be obtained for a special 
case. 


Extremum Char Weight 


When the derivative of the 
char weight with respect to 
time vanishes, the char weight 
is an extremum. It will be seen 
later that the char weight is a 
monotonic function of heat input 
if other conditions are constant. 

Therefore, the extremum may be 
either a maximum or a minimum 
depending on whether the char 
weight is lesser or greater than 
the extremum char weight, as 

shown in sketch (a). Sketch (a) 

From equation (15), the extremum char weight is found to be 



W, 


c.m 


= P 


Ah 


H e " Ml - 


_f Q-r (H e - H v ) (l - f + ari) 


\Ah e ff a Ah eff / 

It should be noted bhat, in general. 


^•c.o 


Ah 


eff 


af 


(17) 


W, 


c.m 


is a function of time. 


The extremum char weight is an important combination of material properties 
and environmental parameters, entirely aside from its interpretation as a lim- 
iting value. The extremum char weight depends on the heating only through the 
ratio of radiative to convective heating rates q r /q c o' and - increases as the 
value of this ratio increases. The volatile fraction and the stream enthalpy 
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strongly influenced the extremum char weight, as shown in figures 2 and 3 . The 
effects of oxygen concentration and radiant heating are discussed later. 



When the char weight approaches 
the extremum char weight, dW c /dt 
approaches zero, from the defini- 
tion of extremum char weight. When 
dW c 

= 0, the rate of pyrolysis is, 
from equation (14-), 


*1 n 


m-n = 


1 - f + aT] H e 


a 


(18) 


w 


For this condition, the rate of 
pyrolysis is independent of k, p, 
e, and Ah e pp. The only material 

property affecting the pyrolysis 
rate is the volatile fraction f. 
This result can be understood when 

dW c 

it is noted that, with — - = 0, 

dt ' 

the rate at which char is formed 
(which is directly proportional to 
the rate of pyrolysis) is equal to 
the rate of char removal. However, 
the rate at which char is removed 
depends only on the rate at which 
oxygen diffuses to the surface. It 
should be noted that the time 

required to approach the extremum char thickness increases as the extremum char 
thickness increases. The extremum char thickness is a function of the material 
properties as well as of the environmental conditions; therefore, the extent to 
which a given test will be influenced by the char weight approaching its extremum 
value is related to the properties of the material. 


Figure 2.- Effect of volatile fraction on extremum 
char weight. 7 \ - O. 75 ; C e = 0.20; J3 = 0.10. 


Caution must be exercised in the interpretation of test results when the 
char-formation characteristics in the test environment do not simulate those 
that exist in the flight environment. When the char thickness approaches its 
extremum value (for example, in low enthalpy tests), the response of the material 
is governed by relations which are quite different from those that apply in 
general. The only material characteristic which affects material performance 
under such conditions is the volatile fraction. 


Solution for Constant Enthalpy 

The most important single heating condition is that for which the stream 
enthalpy is constant and for this condition equation (15) can be integrated in 
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closed form, provided q. r /lc,o an(i c e are constant. All of these conditions 
are typical of tests conducted in ground test facilities. 


Equation ( 15 ) can be simplified 
to the following form by using equa- 
tion ( 17 ): 

a A^c,m W c \ 1 

dWc = 1 + an\l Tj H e - 

dQ c W c t i ~ ri(a Ah + H~ - Hy J 
0 Ah e ff(l + arj) 

(19) 


If the stream enthalpy is constant and 
any radiant heating is proportional to 
the convective heating, equation ( 19 ) 
is a linear first-order differential 
equation in Q c with W c as the inde- 
pendent variable. The solution of this 
equation is 


w c,i - w c 


1 - 


W r 


W, 


c,m 


- 1 + 


Ap 


W, 


c.m 


log. 


c,m 


W, 


c, 1 
c.m 


Qr 


a ^c 

1 + aT) (H e - H W )W C ^! 


(20a) 



Figure 3.- Effect of enthalpy on extremum 
char weight. q r = 0 . 


where 


n(a£h + H e - Hj 
^eff^ 1 + aT l) 


(20b) 


The amount of pyrolysis that has occurred can be determined as a function of the 
heat input and the char thickness by integrating equation (l^) 


IDp = 


w c - w c i + — ^9 

1 + at) H e - H w 


1 - f + 


( 21 ) 


a T) 


f(l + arj) 


The weight of material which has been thermally degraded (that is, pyro- 
lyzed to yield gaseous products or char) is, from equations (l6) and (21), 

a Qc 

w c “ w c i + 

0,1 1 + at) H e - 

W = — + W Cjl 22) 


1 - f + arj 
1 + aq 
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Because of the transcendental character of equation (20), the required weight 
of material cannot be expressed explicitly as a function of the heat input. 
However, equation (20) can be used to determine W c as a function of Q c . 
These results can be substituted into equation (22) to calculate the required 
weight . 


Solution for No Char Removal 


For the case in which no char removal occurs (a = 0), equation (15) has 
the following form: 


dw c 

W c + 1 + T l( H e “ H v ) 

l 

iQ 

O 

V. 

0 

+ 

hi 

dt 

0 Ah e ff 

^eff 


(23) 


This form of the equation would apply to a material having characteristics 
similar to a charring ablator except that the surface does not oxidize and is 
not removed by mechanical forces. An impregnated ceramic shield might provide 
an example of this type of performance. This equation would also apply to 
reentry into a planetary atmosphere containing no oxygen or other chemical 
species which react with the char. 


The solution of equation (23) is 

fw, 


He 

p 


' c -^ + 1 + 


11 ( H e " H wj] + 2(1 - f) ( Q c + Q r) 

f p Ah eff 


.V2 


S 2 I. + 1 + 


Ah e ff J 

l( s a - \) 


Ah 


eff 


(2lki) 


and the total weight of degraded material is 


W = W 


c,i 



(2irt>) 


If an average value of the enthalpy during reentry is assumed, this result is 
identical to that presented in reference 1. 


DISCUSSION 


Certain characteristic features of char formation and removal have been 
noted previously. In general, however, the equations governing the performance 
of charring ablators are too complex to yield significant information from a 
cursory examination. The effects of material properties and environmental 
parameters on system performance can best be evaluated on the basis of a para- 
metric study. All the results presented in figures 4 to 9 are based on heating 
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in an airstream. Therefore, with a given heat input but with different enthal- 
pies, the difference between the two results indicates the extent to which a 
test at one condition fails to simulate performance at the other condition; that 
is, it indicates the difference which results from failure to consider the varia- 
tions in the two test conditions. 


The relation between W/S and — is independent of the value 

' H P Ah ef f 

of p. However, the required weight itself is not independent of p. For 
example, if p and Q both increase by a factor of 2 , the required weight also 
increases by a factor of 2. Little information is available concerning the 
appropriate value of p. Based on a comparison of calculated results and a 
limited number of test results, the value of p for a 5 O- 5 O mixture of phenolic 
and nylon was found to be about 0 . 10 . 


Thermal-Protection Weight 

The conditions for which equation ( 19 ) has been integrated are not repre- 
sentative of an actual spacecraft reentry. During such a reentry, the enthalpy 
decreases from its initial value to a negligible value; however, a significant 
fraction of the total heat input may be experienced before the decrease in 
enthalpy is appreciable. In tests conducted In ground facilities, a constant- 
enthalpy level is normally maintained and the present solution should be valid. 
Therefore, weights calculated from equations (20) and (22) are useful in eval- 
uating the effects of environmental parameters and material properties on 
thermal- protect ion requirements and also in determining the significance of 
ground test results as they apply to a flight environment. 

The thermal-protection weight (eq. (22)) is plotted as a function of the 
total heat Input in figure 4. Although this weight does not include any allow- 
ance for insulation or for factors of safety, it will be referred to as the 
thermal-protection weight. The required weight is higher at the low enthalpy 
level than at the high enthalpy level for all values of total heat input. This 
condition results from the greater effect of oxidation and the lesser effect of 
blocking at the lower enthalpy. 

Over the most important range of heat input, the optimum value of f is 
between 0.4 and 0.6, the lower value of f being associated with the higher 
heat input. However, for values of f between 0.4 and 0.6, the required weight 
is not extremely sensitive to the particular value that is used. 

The effect of total heat input on required weight is shown in figure 5 for 
a value f = 0.4. With this linear total heat input scale, the decrease in the 
slope of the weight curve with increasing total heat input is readily observed. 
With no radiant heat input, the required weight is more strongly influenced by 
the enthalpy level. If one-half of the total heat input is radiant 
(l r /q. Cj> o = l)j the required weight at a constant -enthalpy level is less than 

that for similar total heating in the absence of radiation. This effect is much 
more pronounced at the low-enthalpy level. Available test results, which are 
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Figure 4.- Effect of heat input on required weight of ablating material. 


reported on the basis of temperature rise rather than on the weight of material 
degraded, indicate that the effect of radiant heating on material performance 



is similar to an equal amount of con- 
vective heating. 

The effect of the enthalpy level on 
required weight is shown explicitly in 
figure 6 for selected values of total 
heat input. At the lower enthalpies^ 
the required weight is very sensitive to 
the enthalpy value whereas at higher 
enthalpies the required weight is rela- 
tively insensitive to the enthalpy. 


Thermal-Protection Effectiveness 

It Is difficult to define a simple 
parameter of considerable generality 
which provides a measure of the effec- 
tiveness of char-forming materials in 
the sense that the heat of ablation pro- 
vides such a measure for subliming 


Figure 5 .- Effect of total heating on required 
weight of ablating material, f = 0.4. 


ablators. A parameter consisting 
of a total heat input divided by 
the thermal-protection weight 
required to accommodate this heat 
input, which has been used in a 
number of investigations, appears 
to provide the best available 
characterization of charring 
ablators : 

Q c 

E “ f + Q r 


The weight in this equation is 
that required to assure the back 
surface of a shield subjected to 
the total heat input Q c + Q r 

will not exceed the temperature 
of pyrolysis. Some results of 
the present study are presented 
in terms of the ratio of this 
parameter to ^ n figure 7 . 

The effectiveness increases with increasing stream enthalpy and is greater 
for high total heat input than for lower heat input. The values shown in fig- 
ure 7 are for a value of f of 0.4 which is less than the optimum value at high 
enthalpies. Therefore the estimated system effectiveness in figure 7 is less 
than the maximum theoretical value. At higher enthalpies the effectiveness of 
charring ablators is 8 to 14 times 
greater than the effective heat of 
pyrolysis for the range of impor- 
tant heat loads . 


Accumulation of Char 

The extremum char weight is a 
significant parameter in the per- 
formance of charring ablators; 
however, it is a limiting value 
and not necessarily a value that 
will be attained or even approached 
during a given heating cycle. For 
example, with a low volatile frac- 
tion and a high enthalpy level, 
extremum char weights in excess of 
20 pounds per square foot (assuming 
(3 = 0.10) are predicted from fig- 
ure 2. This value is greater than 
the estimated total thermal 



Figure 7.- Effect of enthalpy level on effective- 
ness of ablating material, f = 0.4; q r =0. 


6 

v, 

lb/sq. ft 



Figure 6.- Effect of enthalpy level on required 
weight of ablating material. f = 0.4; q p = 0. 
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protection weight requirements for the stagnation point of a vehicle entering 
the earth f s atmosphere at hyperbolic velocity. (See ref. 4.) Therefore, it 
is necessary to examine actual char weights in relation to extremum char weight 

as a function of the total 



(a) Low enthalpy. 


- Hg 
^eff 


2 . 



n e ■“ nrr 

(b) High enthalpy, —r = 10. 

^eff 

Figure 8.- Effect of heat input on the relative magni- 
tudes of char weight and extremum char weight. 


heat input . The char weight 
as a fraction of the extremum 
char weight is plotted in 
figure 8 as a function of 
total heat input, with zero 
initial char weight. 

At the lower enthalpy 
level (fig. 8(a)), W c 
approaches W c max at rela- 
tively low total heat inputs, 
particularly when the vola- 
tile fraction is high. At 
the higher enthalpy level 
(fig. 8(b)), the char thick- 
ness is considerably less 
than its maximum value, even 
with high total heat inputs 
and a high volatile fraction. 
Therefore, at low enthalpy 
levels the performance of a 
charring ablator is strongly 
influenced by the char thick- 
ness approaching its maximum 
value, whereas at high 
enthalpy levels, the approach 
of the char weight to a max- 
imum value is barely evident. 
When radiant heating is pres- 
ent, the char weight as a 
fraction of its maximum value 
is greatly reduced for given 
total heat input. 

That W c yAtf c ,max -*- s a 
dominant factor in the per- 
formance of charring ablators 
is evident from the basic 
equation (eq. ( 20 )). The 
total heat input required to 
produce a char weight (with 
W c i = 0), that is, within 
l/e of the maximum char 
weight is plotted in fig- 
ure 9 as a function of f. 

An increase in the enthalpy 
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n 


level by a factor of 5 results in 
an increase by a factor of 20 in 
the heat input required to achieve 
the given value of Wq/W c ^ max . 

LIMITATIONS AND EXTENSIONS 
OF ANALYSIS 


The method of analysis pre- 
sented herein is based on a dif- 
fusion controlled oxidation of the 
char surface. This mechanism has 
been verified in a series of tests 
of 50 -percent -phenolic — 50 -percent- 
nylon specimens. These tests were 
conducted in an atmospheric- 
pressure arc-heated jet with a 
heating rate of about 
l80 Btu/ft ^-sec. For tests of 
other materials or for tests in a 
different environment, other 
mechanisms of char removal may be 
Figure 9»- Heat input required for char weight important. Even when Oxidation 

to reach l/e of extremum char weight. ± s the only mechanism involved in 

W T 

= 1 - q r = °- char removal, the process will be 

c ' m reaction- rate controlled rather 

than diffusion controlled when 
the heating rates, and conse- 
quently the surface temperatures, are sufficiently low. If mechanical forces 
such as shear stresses, differential pressures, or thermal stresses are a fac- 
tor in char removal, these effects must be superposed on the effects of 
oxidation. 

In the present analysis it has been assumed that pyrolysis occurs at the 
desired back surface temperature and therefore no additional insulation or 
cooling is required. Since this condition is not usually satisfied with actual 
materials, the analysis can be used only to calculate the amount of material 
that will be pyrolyzed. The effect of the heat conducted from the pyrolysis 
zone into the uncharred material on the rate of pyrolysis can be accounted for 
by adjusting the heat of pyrolysis to include the heat required to raise the 
material from its initial temperature to the pyrolysis temperature. The actual 
temperature response of the uncharred material and of the structure must be 
determined on the basis of a separate analysis. 

In deriving the basic differential equation, no assumptions have been made 
regarding the functional dependence of the enthalpy or of the radiant heating; 
therefore, if these parameters are expressed as functions of Q c , equation ( 19 ) 
is still a valid relation, and solutions can be obtained by numerical methods. 



IT 


In particular, if the enthalpy and radiant heating are slowly varying func- 
tions of Q c , a solution can be obtained by successive applications of equa- 
tions (20a) and (22); that is, the final value of W c in each calculation 
becomes the initial value in the succeeding calculation. The appropriate cur- 
rent values of al 1 the environmental parameters must be used in each step of 
the calculation. 


CONCLUDING REMARKS 

4 


A diffusion- controlled oxidation mechanism for the removal of char at the 
surface of a char- forming material can he incorporated in an approximate anal- 
ysis, and the performance of such materials in various thermal environments 
can be evaluated- For a given material and environment, there exists an 
extremum char weight. The performance of a material at any time is a function 
of the ratio of the actual char weight to the extremum char weight. 

When the actual char weight approaches the extremum char weight, the rate 
of pyrolysis approaches a limiting value which is independent of the thermal 
conductivity and surface emissivity of the char and of the heat of pyrolysis. 
The extremum char weight increases linearly with the enthalpy level; therefore, 
the actual char weight approaches the extremum value much more quickly in tests 
conducted at low-enthalpy levels than in a high-enthalpy reentry environment. 
This factor must be considered in evaluating the significance of results 
obtained in facilities which provide a low-enthalpy test stream. 

The effects of radiant heating on the performance of charring ablators 
appear in the governing equations only in the extremum char thickness rela- 
tion. The an aly sis indicates that with constant total heat input the effec- 
tiveness of such materials increases with the substitution of a small component 
of radiant heating for convective heating. This increase in effectiveness is 
much more pronounced at low enthalpies than at high enthalpies. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 3 , 1963* 
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APPENDIX 


ANALYSIS OF OXIDATION 


The oxidation of carbon surfaces in hot airstreams has been extensively 
studied analytically. (See refs. 5 to 7-) It is found that at low surface 
temperatures the rate of oxidation of carbon is controlled by the reaction 
rates. At higher surface temperatures, however, the process is diffusion con- 
trolled; that is, the rate of oxidation is limited by the quantity of oxygen 
diffusing to the surface through the boundary layer. For porous carbon sur- 
faces, the temperature at which the process becomes diffusion controlled is 
lower than it is for nonporous surfaces. 

The char layer that is formed when charring ablators are heated tends to 
be porous, as evidenced by the fact that the gaseous products of pyrolysis 
transpire through it. In addition, injection of these pyrolysis products into 
the boundary layer reduces the quantity of oxygen diffusing to the surface. It 
is anticipated, therefore, that the oxidation of the char surface will be dif- 
fusion controlled, even at relatively low temperatures. 

The rate at which oxygen diffuses to the surface is determined from the 
boundary-layer equation for conservation of oxygen. This equation is (ref. 8) 


p’u ^ + p'v = jL/jnag.') 
dx by by \N Sc by) 


The energy equation can be expressed i 


•? •Tl *4* V* ”1 T-l VMT 




(Al) 


p. v gu a 

dx ay dy 


_n_ M + Jifi 

N pr Sy 2 


1 \du 2 

Np r /dy 


(A2) 


Equations (Al) and (A2) are based on a nonreacting mixture of two gases having 
the same heat capacity. If Np r =1 or if bu^/by is small compared with 

b&/by } which is valid at stagnation regions, equation (A2) may be expressed 
approximately as follows: 


P . U 2H + p. T M , |_Lm_ M 

Bx by by \Np r by 


(A3) 


The maximum rate of diffusion of oxygen to the surface is obtained when 
the concentration of oxygen at t 
equations (Al) and (A3) are 


wall vanishes. 

The boundary < 

o 

n 

o 

(y 

H = H w 

(y 


(A4) 
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(A5) 


The oxygen concentration and the enthalpy can be replaced by the dimen- 
sionless variables: 


c 


(A 6 a) 


- H - E 
H = 


■w 


H e “ % 


(A 6 b) 


In terms of these new variables, equations (Al), (A3), (A4), and (A 5 ) are 

1 + -- f - Ifc I) 


. Sh . 1 Sh _ S / (j. 

P U c)x P V Sy 8y \Npr 8yy 


C = H = 0 
C = H -> 1 


(y - 0) 

(y = «) 


(A7b) 

(A7c) 

(A7d) 


Because of the formal similarity of equations (A7a) and (A7b) and their 
boundary conditions (A 7 c) and (A 7 d), if the solution of equation (A 7b) is 


H = 9( x .y^Hp r ) 

the solution of equation (A 7 a) is 

c - T( x ,y,% c ) 


(A 8 a) 


(A 8 b) 


where cp is the same function in both cases. It is shown in reference 10 
that cp(x,y,Np r ) can be expressed as 


cp(x,y,N Pr ) = NprV]_(x,y) 


(a8c) 


Therefore, 


cp(x,y,W Sc ) = N^Vi(x,y) 


(A 8 d) 
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The rate at which oxygen diffuses to the wall is 


& ( 0 2) ‘ H) 


(A9a) 


The aerodynamic heating rate is 




(A9b) 


The ratio of rate of diffusion of oxygen to the surface to the heat-transfer 
rate is 
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